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I discuss the derivation and apphcations of Moving NRQCD (MNRQCD), a generalization of NRQCD which al- 
lows the treatment of heavy quarks moving with a finite velocity. This formalism is vital in reducing discretization 
errors in calculations of large recoil decays, such as _B — > K* + 7 or the Isgur-Wise function at large w. 



1. INTRODUCTION 

Many of the most interesting experimental 
probes of the Standard Model involve exclusive 
weak decays of B mesons. Because the B is much 
heavier than charmed and strange mesons, these 
decay products have large recoil momenta. This 
presents a problem for lattice calculations of these 
processes, since (naively) the spatial lattice spac- 
ing must be fine enough to control discretization 
errors of order (aprecou)- For example, in the pro- 
cess B K* + 7, the K* is kinematically con- 
strained to have a momentum of about 2.5GeV 
relative to the B; this means that an inverse lat- 
tice spacing of lOGeV would be required just to 
reduce (aprecou) to about 25%. 

A closely related problem has already been 
dealt with in lattice studies of B physics, namely 
that the total energy of the h quark is much larger 
than its (dynamically important) recoil momen- 
tum and kinetic energy within the meson. Sev- 
eral solutions have been found; the NRQCD and 
FNAL approaches are the most widely used. 
In this talk I will discuss extending the NRQCD 
method to heavy quarks moving with finite ve- 
locity; this allows much of the recoil momentum 
to be removed from the calculation before dis- 
cretization. The idea is that "removing the rest 
mass" in NRQCD is equivalent to shifting the 4- 
momentum of the b quark, Pf, by an amount pro- 
portional to the time 4- vector, while MNRQCD 
shifts the h quark's 4-momentum by a multiple of 
an arbitrary time-like 4- vector (preferably chosen 



to be the 4- velocity, C/s , of the B meson) 

NRQCD : Ft ^ Pb + rribi, 
MNRQCD : Pb Pb + nibUB, 
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where pt is the shifted 4-momentum of the quark, 
which is relevant for lattice discretization er- 
rors, and nib is an arbitrarily chosen energy shift 
parameter (usually chosen non-relativistically so 
that the kinetic and static masses of the physical 
state are equal) .0 li Ub — i the two methods are 
identical, but when Ub ^ i the MNRQCD shift 
reduces the large spatial components of Pb ■ 

Some work on this method has appeared previ- 
ously. In the same way that MNRQCD is a gen- 
eralization of NRQCD, the "moving static" for- 
malism of Mandula and Ogilvie |^ (MO) general- 
izes the static theory, which in turn is the infinite 
mass limit of NRQCD. Continuum HQET deriva- 
tions are usually done in a Lorentz-covariant man- 
ner; the terms in the continuum limit of the MN- 
RQCD action at some order in 1/M should agree 
(up to field redefinitions) with the corresponding 
HQET action. Finally, Hashimoto and Matsu- 
furu g (HM) have derived the 0{p^/M) terms in 
the MNRQCD action and done numerical stud- 
ies at low recoil; Both MO and HM have also 
studied non-perturbative renormalization of the 
input 4-velocity in the moving static and MN- 
RQCD actions, respectively. In this talk, I dis- 
cuss the kinematics of specific high-recoil decays 
and the expected improvement in discretization 
errors from using MNRQCD. I sketch a method 
of deriving the tree-level MNRQCD action by us- 



^Note that |tp = \Ub\'^ = —1 even when working with the 
Wick-rotated EucHdean theory. 
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ing the Foldy-Wouthuysen-Tani (FWT) transfor- 
mation and use this method to write down aU 
0{\/M) terms in continuum MNRQCD. 

2. DECAY KINEMATICS 

The most ilhistrative decay to consider \s B ^ 
K* + 7. Because it is a two particle decay, the 
recoil 3-momentum of the K* (in the rest frame of 
the B) is fixed to be 2.5GeV, which corresponds 
to an energy of 2.6GeV and a relativistic T fac- 
tor of 3.0. The energy scale fiHAD which controls 
discretization errors in Clover light hadron spec- 
troscopy is a few hundred MeV; inverse lattice 
spacings of 500-lOOOMcV give accuracies of a few 
percent Q. In the B's rest frame, however, 
will be controlled by Pk* , the K* 's 4- momentum. 
This scale is about an order of magnitude larger, 
so inverse lattice spacings of 5-lOGeV might be 
required to achieve few percent accuracies. 

In contrast, consider the rest frame of the 
K* . If we make this the rest frame of our lat- 
tice discretization, then fix* will be the usual 
few hundred MeV, but the B will be boosted 
to r = 3.0; B meson discretization errors will 
dominate. There are two sorts of discretization 
errors we need to consider: "boosted muck" and 
"moving mass" . The boosted muck errors arise 
from Lorentz contraction; the brown muck cloud 
around the b quark will be Lorentz contracted 
in the direction of motion when the B meson 
is moving. This should result in discretization 
errors with a typical scale of TB^J,muck, where 
Hmuck ~ HHAD is the ^ ioT & B meson at rest. 

The moving mass errors are the same ones 
we encountered when boosting the K*] a par- 
ticle of mass TO has spatial momentum |p| — 
TO|v|r = TOvT^ — 1 when moving with velocity 

V and boost factor F. After performing the MN- 
RQCD shift of the 4- momentum in Eq. (0), the 
relevant mass is = — to^, for which most 
NRQCD calculations obtain a value of about 
IGeV. The moving mass discretization scale in 
this case then is A^VF^ — 1 w 3GeV, i.e. roughly 
the same as for the K* in the original frame. This 
is due to the numerical accident that m*j^ « A^; 

V and F don't care which meson is at rest so equal 
"mass" mesons have equal momenta. 



The problem in both the above frames is that 
the discretization scales of the two mesons are 
very different. An optimal frame is one in which 
the two mesons are equally well discretized. Since 
to|^ « As , this means that the two mesons should 
have roughly equal boost factors, i.e. F^ = 
Vk' = yJ{l + Ttot)/2 « 72, where Tb and F^. 
are the meson boosts relative to the lattice frame 
and F(ot, their boost relative to each other, has 
been take to be about 3. This reduction in F 
by a factor of 2 leads to an even bigger improve- 
ment in the spatial momentum, which depends 
upon F like VF^ — 1. The spatial momenta of 
the mesons in this frame are y>k* ~ ps w IGeV, 
to be compared with the previous values of 2.5- 
3GeV. The boosted muck errors in this frame are 
less than twice fin ad ; moving mass errors should 
dominate. Note that this is a worst-case analysis; 
the moving mass errors of the B meson might be 
further reduced by using a different choice of rrib. 

The other example I will consider is _B ^ 
D* + e + h'. Since this is a three particle decay, 
the D*'s recoil momentum can range between 
and 2.3GeV. The high recoil region gives much 
cleaner experimental results for e~^e~ colliders 
running at the T(4S'), while most of the theo- 
retical work has been done near zero recoil. If we 
work in the rest frame of the B, then at maximum 
recoil the D* has 3-momentum, energy, and F of 
2.3GeV, 3.0GeV, and 1.1, respectively. If we use 
MNRQCD, however, we can subtract from the 4- 
momenta of both the B and D* mesons. Since 
A should be the same for the two mesons (i.e. 
about IGeV), again we want to choose the frame 
where they have equal boost factors (which will 
be about 1.025). This implies a moving mass mo- 
mentum of 220MeV, which is about the same size 
as the (not very) boosted muck scale. This means 
that, for a given lattice spacing, these runs should 
have the same discretization errors as light spec- 
troscopy calculations. Even if we decide not to 
use MNRQCD for the charm quark, the Z?*'s spa- 
tial momentum will still only be about 500MeV. 

3. DERIVATION 

NRQCD is a (non-renormalizable) effective 
field theory which reproduces the QCD heavy 
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quark action. As such, one should perform the 
nornial matching procedure to adjust the coeffi- 
cients in the action. When one is working at tree- 
level, however, a (much simpler) classical deriva- 
tion can be used; this is just the FWT trans- 
formation. The FWT transformation consists of 
going to a Dirac basis which diagonalizes 70 and 
then, order by order in p/m, block diagonaliz- 
ing the fermion kernel into non-interacting quark 
and anti-quark sectors. This is equivalent to re- 
quiring that 7o commute with the fermion kernel. 
The 4-momentum is then shifted by rescaling the 
fermion fields by factors of exp (rubjot) . 

To derive the tree-level MNRQCD action, one 
just needs to rewrite the NRQCD derivation in 
covariant language. This is done by writing the 
time components of 4- vectors as their dot product 
with the unit time vector (i.e. 70 is really f), while 
3-vectors are obtained with the "spatial metric" 
s^" = gf"" + f't" (note that h projects onto the 
spatial directions transverse to i, since |ip = —1). 
The p^/2mo term in the NRQCD action, for ex- 
ample, can be written as p^s^^p" /'^rriQ- The MN- 
RQCD derivation then consists of going to a basis 
in which p ^ rather than ^, is diagonal and repeat- 
ing the FWT transformation so that the fermion 
kernel commutes with p . The action obtained 
will be the same as the NRQCD action (written 
in covariant form) with t replaced everywhere by 
U. Since the action kernel commutes with p , 
the moving Pauli spinors can now be rescaled by 
exp^m^pU -x), which effects the subtraction in 
Eq. (^. Dirac spinors (for use in currents) are 
recovered by performing the inverse fermion op- 
erations in reverse order: unrescaling the fields, 
performing an inverse FWT transformation, and 
changing back into a Dirac basis in which 70 is 
diagonal (this last is equivalent to boosting the 
fermion fields). 

There are several subtleties in this proce- 
dure. First of all, the kinetic term, p^{5^^'^ + 
U^U'^)p,y, contains both second temporal and 
spatio-temporal derivatives, which could lead to 
spurious poles in the lattice dispersion relation. 
The solution is a field redefinition 0]; the equa- 
tions of motion are repeatedly used to eliminate 
time derivatives in higher-derivative terms in fa- 



vor of space derivatives. This procedure has been 
carried out in fl; the result is quite nice: 
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where = U'^/Tjy is the velocity of the particle 
and i, j run over spatial directions. This is a 
Lorentz contraction term; when |v| — 1 this term 
projects out the transverse momenta. 

The other subtlety involves the other 0{l/M) 
term in the NRQCD action: a-B/2mo. The tran- 
scribed expression is in terms of ct^^ , the commu- 
tator of the moving frame 7 matrices. The prob- 
lem is that we need to write this term using cTpi/ 
(the commutator of the rest frame 7's) without 
changing the spinor basis. This is achieved by 
boosting only the vector indices. One finds: 



a ■ B 



(3) 



fT • (B + E X v) 



1 + r 



[a . v) (v B) 



where F = F^. The first term is just the moving 
B field written in terms of rest frame E and B, 
while the second term again completes a trans- 
verse projector when |t;| = 1 (the E term is miss- 
ing because v • E x v vanishes). The full 0{l/M) 
continuum MNRQCD action is then 



S 

Hmn 



Q [Dt + Hmn] Q 



D 



2moF 

a, {5^^ -v'v^) (B-l-v X E) 



2mo. 
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